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Abstract 

Using thermal field theory methods, we recalculate axion emission from dense plas- 
mas. We study in particular the Primakoff and the bremsstrahlung processes. The 
Primakoff rate is significantly suppressed at high densities, when the electrons become 
relativistic. However, the bound on the axion-photon coupling, G < 10~^° GeV, is 
unaffected, as it is constrained by the evolution of HB stars, which have low densi- 
ties. In contradistinction, the same relativistic effects enhance the bremsstrahlung 
processes. From the red giants and white dwarfs evolution, we obtain a conservative 
bound on the axion-electron coupling, Qae < 2 x 10"^'^. 
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1 Introduction 



There are recent speculations that axions |T|] might play a role in astrophysics [0, 0. 
Indeed, the cooling rate of the white dwarf G117-B15A has recently been measured 
and is faster than expected from any conventional model. Axions with masses of the 
order of 8(±3) x 10^^/ cos^/3 eV would provide an additional cooling mechanism and 
reconcile models with the data [0. 

On the other hand, it is well known that other constraints on axions come from 
astrophysical systems such as red giants and SN1987A Q]. A rather well-established 
bound comes from the evolution of red giants The upper bound for DFSZ axions 
is rria < 1.1 X 10^^/ cos^/3 eV, quite close to the proposed solution for the cooling 
mechanism of G117-B15A. 

The calculations of the axion emission rates in stellar cores involve a quite compli- 
cated physics (screening and correlation effects). They are based upon kinetic theory 
[0. As a consequence of this complexity, there has been a considerable amount of 
confusion in identifying the correct physics for the photon propagation in dense me- 
dia 1^, P]. The best approach is due to Raffelt |^, |[, and we shall reproduce some of 
his results here. Indeed, given the outstanding implications of the possible existence 
of axions, we propose to recalculate the dominant processes of axion emission in red 
giants and white dwarfs. We shall use a different method, derived from thermal field 



theory un, 11|. The advantage of this method is that it already includes the classical 



limit, but allows also systematic diagrammatic studies of the higher-order effects [12 



The paper is organized as follows: in Sect. 2 we recall some of the basic ingredients 
concerning the spectral density of photons at finite temperature and density. In 
Sect. 3, we investigate the Primakoff process in details. We confirm the previous 
bound on the axion-photon coupling strength. In Sect. 4, we calculate the axion 
emission from electron-nucleus bremsstrahlung through the exchange of transverse 
and longitudinal photons. We discuss our results in Sect. 5. From the red giants 
and white dwarfs evolution, we obtain a conservative bound on the axion-electron 
coupling, gae < 2 X 10"^^. 
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2 The method 



The general framework and methods of gauge field theory at finite temperature and 
density have been discussed extensively elsewhere [Q, and we shall not repeat them 
here. Let us just recall that in the plasma core of dense stars, matter effects are 
properly taken into account by the finite temperature and density corrections to the 
polarization tensor Il{Q). In a heat bath, there exist two modes of propagation of 
electromagnetic waves, transverse (T) and longitudinal (L). In calculating an electro- 
magnetic scattering, it is particularly convenient to make use of the spectral density 

m 

1 

Pt,l(Q) = -Re 



vr g2_nT,L(g) + 2e 
= 5(g2-RenT,L(g))-/?T,L(Q). (1) 

Our convention is to use capital letters only for quadri-momenta, e.g. Q = {uj, q). De- 
pending on whether the photon is time-like or space-like, the above quantity describes 
two different modes: 

i) If the photon is time-like (Q^ > 0), the spectral density turns to a ^-function 
and describes a quasi-particle mode (plasmon). This quasi-particle obeys a dispersion 
relation given, in the non-relativistic limit, by u"^ = uJq + for the transverse mode, 
and by uj'^ = uJq for the longitudinal one. Similar, although much more involved, 
analytic expressions can be derived in the ultra-relativistic limit. Recent works have 
shown that it is possible to avoid the non/ultra-relativistic limit and get closed for- 
mulae The plasmon frequency is given by uJq = e'^Ne/nie for non-degenerate 



and non-relativistic plasmas and by uJq = e Ne/ /i for degenerate plasmas (/i is the rel- 



ativistic electron chemical potential, fi = yPp + 'ml). In the non-relativistic case, the 
longitudinal dispersion relation crosses the light-cone at g = cuq, and eventually con- 
tinues for large values of q until an imaginary part develops in 'n.i^{Q). This peculiar 
feature of a space-like quasi-particle occurs only in a medium. It will be particularly 
relevant for the Primakoff effect, as we shall see in the next section. 

ii) At leading order, the polarization tensor develops an imaginary part only below 
the light-cone, for instance when uj < prq/fJ' in the degenerate case. In the hard-loop 
approximation, the imaginary parts of the transverse and longitudinal components 
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are given by 
ImnT(Q) = 

ImnL(Q) = 
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where the contribution from antiparticles has been dropped. The spectral density 
then turns to some kind of Lorentzian with a complicated structure. One may notice 
that longitudinal and transverse expressions are very similar in the ultra-relativistic 



limit 12 , but not in the non-relativistic case [141. In the non-relativistic and non- 



degenerate case, we recover the standard relation || 
lvaIii^{Q) = ^fKk1)Z e"^^ with 



(3) 



u m 
''"iV2T' 

and is the Debye mass corresponding to the screening of static electric fields. In 
non-degenerate plasmas, = e'^Ne/T, and in degenerate plasmas, kj^ = {e^ /n^)fipF- 
As we shall see in the following, it is important to notice that screening effects are 
absent for transverse modes, while there appears a mass term into the denominator 
of the longitudinal spectral density. More precisely, one has in the static limit 



^2 

RenT(cu^O,g) = 4^ 
and ReUi^iu ^ 0,q) = k^. (4) 
Notice that fc^ also appears in the expression for the transverse propagator. Because 



of the u dependence, one may call this effect a dynamical screening [13]. 

All the above equations have been calculated for the case of electrons or protons. 
Similar expressions can be derived in case of ions, using scalar electrodynamics in- 
stead of quantum electrodynamics. This is legitimate as, in most scatterings, the 
photon does not see the inner structure of the nucleus. Moreover, one does not 
need form factors as we only consider the thermal corrections, which are ultraviolet- 
finite. Equations (3) and (4) are identical, except for the Debye mass, which becomes 
kl = Z^e'^N.onslT. 

The scattering rates are computed either using the cutting rules of Kobes and 
Semenoff ITH] in the real-time formalism []ll|, or by analytic continuation in the 



imaginary-time formalism [T^]. Many examples can be found in the literature, includ- 
ing those discussed here, but mostly in the ultra-relativistic regime [0, |T^, jTU . 



3 The PrimakofF effect 



In a previous work [0, the axion emission through photon exchange was aheady 
studied in detail for the case of an ultra-relativistic plasma. This analysis can be ex- 
tended to the general case of relativistic as well as non-degenerate plasmas. As shown 
in |12[, a systematic study of higher-order effects involves a resummation of a class 



of diagrams. This resummation is best implemented within Thermal Field Theory. 
In this approach, the axion emission rate is calculated through the discontinuity of 
its self-energy. This procedure avoids the complication of using Feynman scattering 
amplitudes, and yields the contribution of several different processes to the axion pro- 
duction. Nevertheless, it remains difficult to compare the results thus obtained with 
kinetic theory calculations 0. 

When we calculate the discontinuity of the axion self energy, three different kinds 
of contributions naturally appear. Out of them, the "pole-cut" term is what we can 
best compare with the classical Primakoff process. It involves the scattering of a 
time-like plasmon with a target (which can either be an ion or an electron) through 
the exchange of a space-like plasmon. To be able to compare our results with previous 
calculations, we shall concentrate on the case where the time-like photon is transverse, 
while the space-like one is longitudinal. The equivalence between Primakoff and pole- 
cut is exact when the photon energies are hard, i.e. larger than the Debye mass or the 
plasmon frequency; it ceases to be valid when the photons are soft, as the space-hke 
photon is affected by Landau damping. 

In our calculation of the axion emission of stars, we shall first compute the dom- 
inant pole-cut contribution, then discuss when the pole-pole contribution becomes 
relevant, and finally find a complete expression valid in the non-degenerate, non- 
relativistic case. 



Following previous work |T2[ , we start from the expression for the pole-cut contri- 
bution to the emission rate of axions 

Tpc{K) = -Y duo smOde qnB{uj){l + nB{uj'))Z2{Q)f3{Q')g{Q,Q'), (5) 

StTK Jujo J-i 

where the same notation as in [0 is used {Q = Q' + K). In the last expression, Q' 
is the momentum of the space-like, longitudinal plasmon, while Q is the momentum 
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of the transverse, on-shell photon (see Fig. 1). For a pseudoscalar couphng one has 
g{Q^ Q') = q^Q'^ sin^ 0, where is the angle between Q and Q' |18| . 



The non-relativistic dispersion relations for a transverse photon give 

cu2 = g2 + ^2 Z2(Q) = 1. (6) 

In fact, the dispersion relation for the transverse photon in the non-relativistic limit 
differs only very slightly from the ultrarelativistic one 0, ^]. We assume that, in the 
calculation of the pole-cut term, the approximation of the dispersion relation given 
above can be used in all regimes to a good accuracy. We need also the 'cut' contribu- 
tion of the spectral density of the longitudinal plasmon. As a first approximation we 
shall take 

where the longitudinal plasmon is static. As the imaginary part of IIl only exists for 
uj' < we expect this approximation to hold well in the case of a non-relativistic 
plasma, where the ion velocity v^, is indeed very small, both in a degenerate and in a 
non-degenerate electron gas. The Debye mass appearing in the photon propagator is 
in fact the sum of the different contributions to the screening (electrons, protons and 
dominantly ions contributions). 

The static limit allows a direct evaluation of 

J dJ{l + nB(o;'))ImnL J ^ImHL ~ nTkl. (8) 

After performing the remaining integration on 6, one gets the final result for the 
pole-cut contribution to the emission rate of axions 



In V - 1 



(9) 



Ak^k^-ujl {^k^ -ujl-kf + kl 

When the axion momentum k is much larger than the plasmon frequency cuq, the 
expression in square brackets simplifies to 



(10) 



The main approximation in the calculation of the pole-cut contribution is the static 
limit for the longitudinal space-like plasmon. This approximation should work well 
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in a large variety of regimes, as it is the interaction with the ions that contributes 
most to both the Debye length and the imaginary part of IIl. The ions being always 
non-relativistic, our approximation is well justified. 

The pole-pole contribution corresponds to the decay of a transverse photon into an 
axion and a plasmon. When calculating it, we must use the longitudinal plasmon dis- 
persion relation. In the low-momentum limit, this dispersion relation gives a constant 
energy ujq, while in the large momentum limit, it tends asymptotically to u = 
Moreover, the Jacobian includes a factor u — f *g, which suppresses the process once 
the asymptote is reached. In the case of a very dense plasma, the asymptote coincides 
with the light cone, and only g's of the order of or smaller than the frequency of the 
plasma contribute. Nevertheless, the phase space increases considerably in the case of 
non-relativistic, dilute plasmas, where q can be as large as tu/f*. Notice that in this 
case we have a space-like solution for the longitudinal particle-like propagation in the 
plasma. 

It is hard to find a good analytical approximation for the pole-pole contribution 
in regions other than non-relativistic and non-degenerate. Therefore, we perform the 



calculation numerically. We start again from the results of |12|, 

rpp{K) = —- / dujnB{uj){l + nB{uj'))Z2{Q')Z,{Q)g{Q,Q'). (11) 

Here, u^q) and Z2{Q) are the same as in Eq. (6), while the longitudinal dispersion 
relations give ^ 



1 UJ ^ oj + Vpq 
m 

2 Vpq' uj' — Vpq' 



and Z2{Q') 



3g'2 \ujl-uj'^ + vlq'^ 



(12) 



We show the result of our numerical calculation in Fig. 2. The pole-pole contribution 
gets large at low densities, and is indeed dominated by the exchange of space-like lon- 
gitudinal plasmons. The curve follows remarkably the pole-cut term, that is Eq. (9). 
At high densities, though, the pole-pole is dominated by the exchange of time-like lon- 
gitudinal plasmons, and is still ~ 50% of the pole-cut rate. However, in this region, 
the rate is dominated by the pole-cut term with exchange of a space-like transverse 
plasmons (we shall comment more in detail on that point below). 
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We have also plotted the production rate predicted by Raffelt. We see that our 
prediction only agrees in the region of low densities. We shall argue that this is 
consistent with the approximations used in both calculations. 

In his calculation, Raffelt used a form factor different from the standard Debye- 
Hiickel one. His justification of this point relies in the stationarity of the targets. 
The relevant time scale for a Primakoff-type event to occur is the time it takes a 
photon to cross a region of the size of the Debye length. If the incident photon is 
very energetic this time will be k'j^ ^ much smaller than v^k]^^ the time an electron 
takes to cross the same region. But whenever the temperature is of the order of cuq, 
the velocity of the average plasmon excitation will be much smaller than the speed 
of light, and the argument above will cease to be correct. Indeed, we see from our 
calculation that when the densities grow, the pole-cut term becomes dominant. In 
fact, it is easily proved that our pole-cut contribution is equivalent to the result that 
would be obtained using the standard Debye-Hiickel form factor in the calculations 
of i- 

Within the non-relativistic, non-degenerate approximation, a compact result for 
the axion emission rate can be deduced, which reproduces exactly the result of ref. 0. 
If we consider both pole-pole and pole-cut contributions together, the axion emission 
rate becomes 

/-oo /•+! 

T{K) = -- duo sinOde qnB{uj){l + nB{uj'))Z2{Q)pM)g{Q^Q')^ (13) 

where pl is now the complete spectral density for the longitudinal plasmon. In terms 
of the electric permittivity e^, one has P 

/?l(Q') = -T^Im- (14) 
In the classical limit, Eq. (13) can be approximated to 

T{K) = —k'^HBik) / duo'— / rf(cos0)sin2 -Im— . (15) 
Stt J-oo uo' J-i \ eL/ 

At this stage, it is quite convenient to use the Kramers- Kronig sum rule 



We are left with the angular integration, which can be solved exactly, leading to the 
final result: 



G_2 
IGtt 



T{K) = —TkluBik) 



(17) 



which indeed coincides with ref. . We notice that in the large-/c region, the sum rules 
are fulfilled by the pole-cut term, while in the \ow-k region the pole-pole dominates. 

Let us now comment on two other points: 
i) The main approximation, which leads to our final result, is the static limit in the 
space-like photon propagator. It is obviously an overestimate. However, it is certainly 



a very good approximation in the non-relativistic limit as < qp/E = gy3T/M ~ 
0.003g (for ^He). As axions mostly originates from ion scatterings, we see that this 
approximation is well justified. 

ii) In Fig. 2, we also show the contribution coming from the exchange of transverse 
photons. Unfortunately, we have not been able to compute analytically the axion 
emission rate in this case, except in a few limiting cases. However, we can use the 



numerical codes developed previously for the ultra- relativistic regime [|18| , |T^. With 



the help of the results of [n, 15 1, the implementation of the general case of a degenerate 
gas of electrons is straightforward. Notice that, in this case, the dominant process 
comes from scatterings where the targets are not ions but electrons. Still, the Debye 
mass that one has to use in the photon propagator (see Eq. (4)) is the sum of all the 
contributions. 

The justification for neglecting transverse space-like excitations is because of their 



usual small fiuctuations in the non-relativistic limit |T^. However, and contrary to 
the longitudinal case, these excitations are not screened and can therefore contribute 
significantly when the Debye mass gets large. Also, as the degeneracy of the electron 
gas increases, the spectral density of transverse excitations becomes as large as the 
longitudinal ones. And this is definitely the case in red giants, when the density 
reaches p = 10^ g/cm^ [|l^. 

The numerical results shown in Fig. 2 clearly demonstrate this feature. 
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4 Bremsstrahlung 



Another important contribution to the emission rate of axions can be the bremsstrah- 
lung process. In this plasmon is absorbed by an electron or ion, and an 
axion is emitted. In contradistinction to a Primakoff-type process, bremsstrahlung 
occurs via the fermion-axion coupling, —ig^'-f^'^^a- A derivative axion coupling was 
introduced by Raffelt, to set right the potentially incorrect consequences of the usual 
pseudoscalar coupling. Nevertheless, in the case of the bremsstrahlung process, it can 
easily be shown that both couplings can be used, and lead to identical results. 

In the case of ultra-degenerate plasmas of electrons, several approximations can be 
performed. On the one hand, Pauli-blocking effects imply that only electrons near the 
Fermi sphere contribute [0. These electrons are very energetic, with energies of the 
order of /i, and the hard-loop approximation will be justified. Another consequence 
of Pauli-blocking is the strong suppression of the e~e~ bremsstrahlung processes in 
comparison to e~-ion ones. 

We shall be concerned with the imaginary part of the polarization tensor, taking 
only into account the exchange with the nuclei in the photon propagator. We only 
consider the space-like photon, as the time-like photon is related with the Compton 
process, which has been shown to be subleading . In a strongly degenerate plasma, 
this is easy to understand, as the plasmon density is automatically damped by a factor 
g-a;o/T_ Using the cutting rules of Kobes and Semenoff, the emission rate is written 
as (see Fig. 3) 

X {e{qo - po) - riFipo - go)) 6 ((P - Q)' - ml) 

X {e{-qo) + nsigo)) f3{Q) E \-^\^ (18) 

spins 

where the amplitude is given by 

A= el{Q)u^\P-K) 

After summing over all spins and photon polarizations, we can split the total squared 
matrix element into its longitudinal and transverse parts. Using the hard-loop ap- 



P2 — rnl 



r 



+ me 



{P-Q-KY- ml 



u'iP-Q). (19) 
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proximation, we obtain, 

E 1-4 It = ^^p^[2(X-gr + gV(-l + (k.qf)]; (20) 

spins V ' 

spins \ / 

In a very dense plasma, Fermi distribution functions can be well approximated by 
step functions. This approximation yields the following constraints on the energies 

II + ko < po < II + qo, (21) 

and allows us to write 

nB{qo)nF{po - go)(l - ?^f(Po - ko)) ~ ns(go)(l + ^5(50 - /^o))(?o - ko)5{po-n). (22) 
At this stage, the integrations over po and p are easily solved, leading to 

^tAK) = ^-fPFj^J^PT,L{Q)nB{qo){l + ns{qo-ko)) 



x{qo - ko)6 ((P - Q)' - ml) E 1^ \t,l- (23) 

spins 

We shall first calculate the transverse contribution to the emission rate, and pro- 
ceed later to the longitudinal case. Before performing the angular integrations, let us 
have a look at the spectral densities of the exchanged photons. The exchange with 
the nuclei involves the space-like contribution of the spectral density (Eq. (1)). In the 
case of heavy, non-relativistic ions we can use the Boltzman statistical factor in the 
calculation of Imllx, and we obtain: 

/tT T ^0 3 

-klt',-e^ni\qo)e~^^* , (24) 

where — 3T/M ^ 1 is the average speed of ions and M is the mass of the nuclei. 
The exponential factor allows us to approximate Qq <S q. Moreover, if we recall 
that screening effects are absent for transverse modes, we see that RcIIt <^ q^ and 
Imllx <S q^. Therefore the transverse spectral density has the following form 

. -ilEHiM, (25, 

TT 
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which just corresponds to using bare propagators. Performing the angular integrations 
is not a trivial task, since three combined angles appear, namely k • q, k • p and p • q. 
To ease the calculation, a usual trick, based on the isotropy of the emission rate, can 
be used, 

T{K) = d^k'TiK'). (26) 

Notice that if k • p for fixed p is not constrained, q • p for fixed q is on the other hand 
bounded. 

We arrive at an expression for F where only one integral is left unsolved, 

9l klvlT 



1927r4 kpl 

X / dqo {l+nB{qo)) 
Jo 



Pf fi-PF \mi ^ 



QQ k 
g 2T g 2T . 



, ^ ^ . (27) 
Finally, we can consider the energy loss rate per volume, related to the emission rate 

by 

1 

Qt^t^ dkk'V^iK). (28) 

Prom the last expression of the emission rate, we can deduce the value of the energy 
loss rate for a transverse photon exchanged. 



9le klvlT^ 



48v/27r6 pI 
/7r^_ 5 
4 3 



Ain^i±^ + 2f4 
Pf I^-Pf \mi 



X 



1 /■°° \ 
C(3) + dx x\l + nB{x)) U{x/2) j . (29) 



We now turn to the longitudinal case. The imaginary part of the polarization 
tensor for heavy, non-relativistic ions is. 



ImnL(Q) = — - n-B\qo)e-^^- 



(30) 



Its structure is such that it allows us to approximate go <S Q, as in the transverse 
case. Then, RellL — kjj, where is the Debye mass. Moreover, ImllL is negligible 
compared with (Q^ — RcnL(Q))^ in the denominator of the spectral density, and the 
spectral density can be written as 

1 ImnL(Q) 



n (g2 + kl) 



2 \2- 



(31) 
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The calculation of the emission rate and energy loss rate, although more involved, 
proceeds in the same manner as in the transverse case. The final result is 



9le 



13^ klT 



327r5V 2 kp^v^ 



Pf 



2 + 1 -^In 



X 



dqa (1 + nsiqQ - ko)) (go 



Pf l^-PF 



ko)e 2T 



q^dq 



2 

% 3 



(g2 + kj 



^.(32) 



'90 -r n.^j 

The calculation of the energy loss rate may be easier, using the Parseval theorem 
between functions and their Fourier transforms: 



dq 



q^ 



go 



(g2 + klY 



2 

_% 3 



Ako Jo 



dxe 



3 ^10 

2 kjjv^ 



1 



/3 XQo 



2 knv^ 



(33) 



Performing first the integration over qQ, and then that over k, an exact expression is 
obtained for the energy loss rate, with the remaining simple integration: 



Qi 



1287r6 



pI 



2M^ + 1 -^In 



mt 



X 



dxe ^ r^F 



\x 



) Pf fJ'-PF 
(C(2,F(x) + l) + F 



+ 3xr, (C(3,F(a;) + l) + F- 




- 9 

- + -xr^F 



-1, 



X] 



(34) 



withF(x) = i(l + y6r. 



and the generalized ^-function ({a,b). 



The longitudinal energy loss rate is obviously much larger than the transverse 
one, by a factor l/v"^. The result can be decomposed into an electronic part and a 
remaining integral, which contains the physics of Landau damping with the nuclei. 
This integral depends on a single parameter, kov^/T, which can vary between 1 and 
20 for the physical cases of interest. Numerically, we find that this integral is about 
11., and is almost independent of r*, within 6% accuracy, so that we might as well 
take its value ioi k^ = {]. This seems to be in agreement with previous works, where 
it was found that the Debye screening is negligible 0. The temperature dependence 
of the rate is therefore as T^. 



5 Discussion 

Using thermal field theory methods, we have calculated the axion emission rate 
through the Primakoff mechanism. We have found that relativistic effects strongly 
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suppress the rate at high densities, compared with the non-relativistic approximation. 
Nevertheless, the axion-photon couphng constraint, G < 10^^° GeV^^ [|, |^, is not 
affected by these results. Indeed, this constraint comes from the evolution of HB 
stars, where the density is about p ~ 10^ g/cm^, and where the electrons are clearly 
non-relativistic. For this type of coupling, axion emission has no effect in RG stars. 

We have also calculated the axion emission through the electron-nucleus bremsstrah- 
lung process, for the case of relativistic, ultra-degenerate electrons. Our final result 
is shown in Eq. (34). For the red giant case, the density is p = 10^ g/cm'^, so that we 
have 

6, = <?i 3.8 X 10^6 f-^) erg/g/s. (35) 



.lO^K. 

The difference between the non-relativistic result and the relativistic one is about a 
factor 2 in favor of the relativistic case. At T = 10*^ K, our result is almost the same 



as the one used in \21 , where a numerical code was used to set the constraint on the 



axion coupling. Hence, our result translates into the following bound 

Qae < 1.5 X 10-^\ (36) 

Should we use the usual limit for the extra energy loss rate that is tolerable, < 
30 erg/g/s |Q, we would get a less restrictive bound, (^ae < 3 x 10^^"^. 

For the axion luminosity in white dwarfs and using the same parameters as in [Q, 
we find 

i„ . ,2. T X 10-Le|; (^)\ (37) 

a result with the same temperature dependence as in previous works, but somewhat 
larger |^ Again, the origin of the discrepancy lies in the relativistic effects. 

Requiring that the axion luminosity should not exceed the photon luminosity, we 
obtain the following bound 

gae < 1.8 X 10-l^ (38) 

which is sligthly weaker than the RG bound. 
In the DFSZ model, Eq. (RBI) gives 



rria < 5.3 x 10" W/ cos^ /3, (39) 
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which is as good as the SN1987A constraint, with a much lower uncertainty. In 
conclusion, we see that the evolution of red giants and white dwarfs leads to roughly 
the same constraint on the axion-electron coupling strength. These results can be 
summarized in a very conservative bound 

<7ae < 2 X 10-'\ (40) 



We conclude on the cooling of G117-B15A. The authors of |0] have used a result |^ 



that apparently underestimates the axion emission rate by a factor 4. Consequently, 
the value of Qae needed to account for the fast cooling of G117-B15A is weaker by a 
factor of 2 compared than their estimate. A value as low as Qae = 1-2 x 10^^'^ would 
reconcile the models with the data. This may have no consequence either on the red 
giants evolution or on SN1987A. 
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Figure Captions 



Fig. 1 Feynman diagram for the Primakoff process. 

Fig. 2 Axion energy loss rate due to the Primakoff effect for a wide range of densities 
typical of red giant stars. The solid curve is our result (Eq. (11)). Also plotted 
are Raffclt's result (dashed curve) and the pole-cut rate with transverse photons 
only (dotted curve). The value of the axion- photon coupling has been fixed to 
G = 10-1° GeV-^ 

Fig. 3 The two cut graphs that arc needed to evaluate the bremsstrahlung contribu- 
tion to the axion emission rate. 
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